Periodic volume fluctuations with infinite horizon: 
Intermittency enhanced Fermi acceleration 
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A particle diffuses in a lattice with periodically forced volume. In the finite horizon case with 
bounded distance between collisions and strongly chaotic dynamics, the velocity growth (Fermi 
acceleration) is well described by a linear Boltzmann equation, leading to an asymptotic universal 
non-Maxwellian velocity distribution scaling as v ~ t. The infinite horizon has intermittent dynamics 
which enhances the acceleration, leading to v ~ tint and a non-universal distribution. 



Periodically forced thermally isolated systems exhibit 
many interesting phenomena from stabilization [JJ to ex- 
ponential acceleration [5] . They are of intense interest for 
trapped atom [3J and ion [¥J experiments, as well as astro- 
physical problems such as transport of comets [3] . Evolv- 
ing in a spatial coordinate, they also describe transport 
in periodic structures [6j. Typically there is unbounded 
growth of the energy, the phenomenon of Fermi accelera- 
tion [7] (FA). Very recently, many researchers have sought 
analytical descriptions of energy distributions in such sys- 
tems [8HTT]. In rather general circumstances a Fokker- 
Planck (FP) equation can be derived, incorporating the 
average and variance of the work per period [5] . The first 
example in [8], extensively investigated elsewhere [5HT5] 
consists of a particle moving freely in a container ( "finite 
billiard") or amongst obstacles ("Lorentz gas") with os- 
cillating boundary but fixed volume. Here we relax the 
fixed volume condition and tame the resulting wild fluc- 
tuations in the energy distribution. We consider a variety 
of regimes distinguished by free paths without collisions 
(infinite horizon, IH) or lack thereof (finite horizon, FH). 
Intermittent IH-type motion is found to actually enhance 
the rate of acceleration, including for the recently intro- 
duced "dynamically infinite horizon" [5], that is, IH for 
only part of each cycle. 

Periodically oscillating billiard(-like) models exhibit- 
ing Fermi acceleration include the ID bouncer [16] and 
stochastic simplified Fermi-Ulam [TU] models. In the lat- 
ter (and often elsewhere), the simplifying assumption of 
the static wall approximation (SWA) was used, where the 
boundaries are fixed (hence trivially having fixed volume) 
but the particle changes its velocity as if it were moving. 
Many oscillating two-dimensional billiards have also been 
considered and lead to Fermi acceleration. It is conjec- 
tured that this includes all chaotic geometries [TU [JS] , as 
well as the ellipse [T7]. The breathing case (fixed shape) 
has been studied in detail [18], leading to slower growth 
of velocity than other typical models. Fermi accelera- 
tion is normally prevented by dissipation in the dynam- 




FIG. 1. Part of the triangular Lorentz gas for infinite horizon 
(IH, left), finite horizon (FH, center), and confined (C, right). 



ics, although scaling laws relating the final energy to the 
strength of the dissipation and other system parameters 
can be observed [IS] . 

Jarzynski and Swiatecki [12) . showed using moments 
that for fixed volume time-dependent billiards, the even- 
tual distribution of velocities is exponential, in con- 
trast to the Gaussian distribution of an equilibrium gas. 
Jarzynski [20] then described an FP equation approach 
for a slowly varying billiard (or fast particle) giving an 
explicit calculation of the rates of increase of the en- 
ergy and its variance. Bouchet, Cecconi and Vulpi- 
ani |21j in an astrophysical context applied a linear Boltz- 
mann equation to obtain an exponential velocity distri- 
bution. More recent innovations have included a hop- 
ping wall approximation replacing the SWA [13;, and a 
Chapman-Kolmogorov equation replacing an FP equa- 
tion [TU]. Here we retain the simpler FP approach, but 
treat the wall collisions exactly. Many of these tech- 
niques are also relevant to stochastically moving bound- 
aries, such as the finite temperature Lorentz gas [TTll2"2"] . 

Our model is a two dimensional Lorentz gas (Fig. [TJ , 
which consists of a collection of circular scatterers in an 
extended domain. Fixed random [53] and periodic [23] 
scatterer arrangements have been widely studied for the 
last century. In the periodic case, FH leads to normal dif- 
fusion (displacements scale as y/i), while IH leads to log- 
arithmic superdiffusion (\/t hit). The collision (discrete 
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time) dynamics is strongly chaotic, with exponential de- 
cay of correlations in both FH and IH cases, however the 
long flights of the IH case leads to C jt decay of correla- 
tions in continuous time. 

We place the scatterers on a triangular lattice, with 
each unit cell having unit area, so the distance between 
the centres of neighbouring scatterers is (4/3) 1 / 4 . The 
triangular Lorentz gas is IH for r < rji = (3/64) 1 / 4 rj 
0.465, FH for r H < r < 77 = (1/12) 1 / 4 0.537, at 
which the scatterers start to intersect, and confined (C) 
for 77 < r < r B = (4/27) 1/4 w 0.620 at which point 
the dynamics is blocked as there is no space outside the 
scatterers. The area available to the billiard particle is 



A(r) = 



1 — irr 2 
6 arccos — 



r < r 1 
r > rj 
(1) 

Here we consider time-dependent scatterers, with ra- 
dius r(t) = R + Asint and velocity u(t) — r'(t) = A cost; 
with suitable scaling this system is equivalent to any 
other with different lattice spacing and frequency. There 
are several scenarios depending on R± = R± A: 



IH: Infinite horizon 

IFH: Infinite, finite horizon 

IFC: Infinite, finite, confined 

FH: Finite horizon 

FC: Finite, confined 

C: Confined 



R+ < r H 
R- < r H < R + < 77 
R- < r H < 77 < R + 
rii < R- < R+ < 77 
th < R- <ri < R+ 
7/ < i?_ 



For the numerical simulations we choose A = 0.03, 
which allows all the above cases except IFC. A square 
lattice Lorentz gas has no finite horizon, and so a time- 
dependent radius would allow IH, C and a further case 
IC allowing infinite horizon and confined but not finite. 

We first discuss Fermi acceleration for the finite or con- 
fined geometries. The billiard particle moves freely, col- 
liding with the scatterers according to [18] 



2n(u — n • v 



(2) 



where v + (v_ ) is the velocity immediately after (before) 
the collision, n is an outward unit vector normal to the 
disk at the point of collision, and we use v± — |v±|. The 
incoming angle 8 with respect to the normal is defined by 
— n • v = cosd. If the particle with V- < u is overtaken 
by the scatterer then 8 > tt/2. We define 8 > so there 
is a 1:1 relation between 9 and the outgoing speed v + ; 
later this requires an extra factor of two. Eq. ([2| gives 

v\ = v + • v + = v 2 _ + Auv- cos 9 + in 2 (3) 

Thus the change in velocity lies between zero and 2u. 
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FIG. 2. The particle moves in the negative x direction. The 
(almost) parallelogram ABDC denotes the initial positions 
that will collide in time St, taking the moving boundary into 
account. 



This system exhibits Fermi acceleration, so almost all 
initial conditions to lead to unbounded velocity; in this 
regime ti> 1 and d»u. Thus the particle is effectively 
in a Lorentz gas with slowly varying radius, and as in the 
static case, exponential decay of correlations. The only 
quantity not randomised by the dynamics at short times 
is v, which is a constant of motion for the static case. 

Thus we may describe the system by a distribution 
function f(v,i)5v giving the probability of observing the 
particle with speed in the interval [v, v + Sv] at time t, 
hence normalised so J Q f(v, t)dv — 1 for all t. The prob- 
ability of finding the particle in a region of the full phase 
is, under this assumption, f(v, t)5v ^ ^£^L w here ip de- 
notes the direction of the velocity, including the relevant 
normalisation factors. Here, and often later, the time 
dependence of r (and hence A) has been suppressed. 

The distribution f(v,t) evolves due to collisions with 
the scatterers, which make small changes of order u to the 
speed. The collisions depend on one distribution function 
and the known position and velocity of the scatterers, 
so the treatment here is a linear Boltzmann equation, 
similar to Ref. [21 . This equation ignores correlations 
between collisions, however we will include their contri- 
butions later using the results of Ref. [20] . 

The general form of a linear Boltzmann equation is 

ft = J \p(v, v')f(v', t) - p(v', v)f(v, t)]dv' (4) 

where subscript t (and later v) is the partial derivative 
and p(v,v') gives the probability density for a collision 
taking v' to v. In our case, we need to compute the 
probability of a collision taking v_ e [v' ' ,v' + Sv 1 ] to 
v + G [v,v + Sv] at a time in [t,t + St] by integrating 
the distribution over the set of trajectories with the ap- 
propriate collision. 

For r < rj the cross section is independent of ip, 
so we take ip — 7r; larger radii or different chaotic 
maps would need to take the ^"dependence into ac- 
count. The trajectory hits the scatterer at time t at a 
point A(r(t) cos(9, r(t) sin#); See Fig. [2] A trajectory 
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hitting the scatterer at angle 9 + 59 at time t reaches 
it at B(r(t)cos(9 + 59),r(t)sm(9 + 89). To reach the 
scatterer at t + St, the particle at time t will be at 
C(r(t+5t) cos 0+V-6t, r(t+5t) sin0) or D{r{t+6t) cos(#+ 
59) + vSt, r(t + St) sin(0 + 59)) respectively. 

We need only leading order in the perturbations, 
so ABDC is a parallelogram, with area r(t)(u + 
t>_ cos 9)595t. We integrate over tp but divide by its nor- 
malising 2-7T. Thus 



P(v+,v-) 



2r , n ,f dd \ 

■ — (v- cos 9 + u) \ - — 



(5) 



where the final derivative comes because we used 9 to 
denote the collision variable rather than v + ; they are 
related by Eq. ^ . The factor of two comes from consid- 
ering both positive and negative angles (see above), and 
the minus sign from the sign of the partial derivative. 
Substituting for 9, we find 



p(v + ,V- 



rv + (v? — v_) 



Au^8u 2 (v 2 + +v 2 _)-(v 2 + - v 2 _) 2 - 16u 4 



(6) 

Anticipating the expansion in powers of u, we now write 
v + = v— + 2su so that s ranges in the fixed interval [0, 1], 
and use this in the linear Boltzmann equation 

J^/O, t ) = j i J Q fr(«> v - 2su )f( v - 2su > t) 



—p(v + 2su, v)f(v, t)] 2uds 



(7) 



For very small velocities (v < 2u) we should modify the 
limits of integration to ensure that the arguments of p 
are both positive, however in practice this is not impor- 
tant as we are interested in long times after which the 
distribution is almost all at large velocities. 

We now consider times of order unity, that is, the pe- 
riod of the oscillations. The Boltzmann equation as it 
stands is not tractable, being explicitly time-dependent. 
Noting again that for typical particle velocities v u, 
we expand the right hand side of the Boltzmann equation 
in a power series in u, a Kramers- Moyal expansion |25j . 
The functions / and p are expanded in powers of u, which 
then allows the integral to be performed, leading to 



fi 



r 
A 



-lTu(f + Vf v ) + 



8u 2 vf vv 



0(u 3 



(8) 



This is now used to determine f(v,t) at long times. We 
note that when r < rj, ie in the IH, IF and FH cases, 
A/2 = —txtu is the term that appears in front of the 
first two terms on the right hand side. Presumably this 
term is also Aj2 for r > rj, as in Ref. [20], which gives 
a comparable equation. [26] Note that terms involving u 3 
and higher are significant only for velocities of order \ft 
or less, thus they do not contribute to the main scaling, 
which is order t. 




FIG. 3. Convergence to the distribution, Eq. ( 13 1 for R = 0.5 
(finite horizon). 



We now come to the main issue as regards the fluctu- 
ating volume. During each oscillation, the particle makes 
of order v collisions with the scatterer during each of the 
expanding and contracting phases, thus increasing and 
decreasing its speed by amounts of order v. Based on an 
approximate averaged equation|27| we propose the fol- 
lowing ansatz: 



f(v,t) = a(tf^e-<W t 



(9) 



where a(t) is bounded and 27r-periodic. Substituting into 
Eq. {8} gives an ODE involving the oscillatory r and u: 



. f , nur\ 8 u 2 r , 

a = (t r- a -a 2 

V A J MA 

This is a Bernoulli equation, with solution 



i(t) 



. f* u 2 r , . 
i(t) = - I -^dt' 



A 



Long time behaviour is characterized by 

m _ 4 



I — lim 

t— >oo t 



3tt 



2 " U 2 T 



VA 



dt. 



(10) 



(11) 



(12) 



an elliptic integral depending on R and A that is easy to 
evaluate numerically |28| . Thus our final expression for 
the velocity distribution function is 



nv,t) = ^ 2 e 



-v/(it) 



(13) 



where V = v\[A and F(V, t)dV = f(v, t)dv. More gener- 
ally, in an adiabatic compression we expect the entropy to 
be roughly invariant, and indeed it is just a function of V 
for this two dimensional ideal gas. Thus in thermally iso- 
lated systems with fluctuating volume the entropy, rather 
than energy or velocity, is the primary variable. 

This distribution is confirmed numerically in Fig. [3] for 
R = 0.5 and A = 0.03, so the FH regime. For this case 
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numerical integration gives I = 0.001327658. We simu- 
late the full (not SWA) time-periodic Lorentz gas. De- 
termining the time until the next collision thus involves 
the solution of a transcendental equation, using the ro- 
bust quadratically convergent method proposed in [IB] , 
A million initial conditions of particles are chosen using 
a Maxwell-Boltzmann distribution at a temperature of 
10~ 4 consistent with \u\ < A = 0.03. Particles which are 
very slow may not collide during the simulation time, and 
so delay convergence to the limiting distribution. Note 
that A more than doubles over the cycle and would be 
clearly visible in Fig. [3] if not incorporated correctly. 

Next, we consider infinite horizon. Here, the correla- 
tions are non-integrable, so the Boltzmann equation can- 
not be applied, and in general the velocity distribution 
is non-universal. We will however determine the scal- 
ing of velocity with time. There are two cases, IH (pure 
infinite horizon) and IF (infinite-finite, also called dy- 
namically infinite). In the IF case the time of free flights 
is bounded by the period, but since the velocity can be 
arbitrarily large, the distance is not bounded. In spite 
of the long flights, most particles at any given time have 
collided much more recently than an oscillation period, so 
the velocity distribution is still approximately a function 
of V as before (see Eqs. [I|[l3|). 

For both IH and IF we expect the Fermi acceleration 
to be of the form v ~ tint, as follows: Consider a flight 
of time t between successive collisions. The probability a 
flight of length between r and r + dr is (neglecting mul- 
tiplicative constants) of order (vt)~ 3 cI(vt) — v~ 2 T~ 3 dr 
for r ^> the typical flight time. See eg Ref. [M] for 
computation of the constant in the static case. 

Collisions normally occur with a rate f=a v, so a long 
flight avoids a change of velocity « vt to the particle, 
noting that u is at most of order unity. Thus the pertur- 
bation 8 In v = Sv/v is of order t. The effects are however 
effectively uncorrelated, so we must add variances, which 
appear in proportion to their probability 



(Slav) 



-\ 2 dr 



(14) 



The lower limit of integration is the typical time u" 1 and 
the upper can be taken as the largest time found in the 
trajectory, however there is no extra velocity perturba- 
tion for free paths of time greater than the period (of 
order 1). Thus we find that per collision the variance 
scales as lnu/w 2 . The number of collisions required to 
reach paths of order unity is about v 2 / Inv, which is less 
than the simulation time, noting that in the finite horizon 
case the Fermi acceleration is of order v <~ t. 

The particle thus undergoes a random walk in In?;, 
taking a number of collisions v 2 /lnt>, and hence a time 
v/ In v to take each step. The total time for the trajectory 
is dominated by the largest value of v in the path, so that 
the velocity is typically of order t Int. This argument fol- 
lows through for both infinite (I) and infinite/finite (IF) 
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R=0.12 
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FIG. 4. Velocity growth for different R. The curves start 
horizonal (v ~ t), and those with R < 0.45, where the infinite 
horizon is noticeable, then increase linearly (v ~ tint). The 
lowest curve, R = 0.58, is a billiard completely confined for 
all time, exhibiting the weakest Fermi acceleration. 



cases, although it is more pronounced (larger coefficient) 
in the former. Note also the paradoxical effect in which 
the intermittency leading to a long time without colli- 
sions is actually responsible for increasing the (collision- 
driven) Fermi acceleration. 

The dependence of the Fermi acceleration on the radius 
(and hence the finite/infinite horizon status) is shown 
in Fig. [4j Here the quantity plotted is the root mean 
square velocity. In time, the regimes are (a) dominated 
by the initial Maxwell-Boltzmann distribution, (b) lin- 
ear growth of v, (c) for IH, increase as the logarithmic 
Eq. ( 14 ) starts to dominate the normal linear accelera- 
tion. 



The transitions depend on the coefficients of the 
relevant terms, moving to shorter times as the radius is 
decreased. 

The FC parameter range, in which the particle is al- 
ternately confined and unconfined, needs further study. 
In particular, it is interesting to investigate the many 
rapid collisions undertaken by a particle near where the 
two scatterers touch, a new "dynamical cusp" mode of 
intermittency. The possibility of an unbounded number 
of collisions suggests that since in each (approximately 
perpendicular) collision, a fixed quantity u is added to 
the velocity, the velocity itself can become unbounded in 
finite time for a small set of initial conditions, a further 
example of intermittency enhanced acceleration. 
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